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Abstract 

We study neutral and charged Higgs boson production in association with stop and 
sbottom squarks at the Large Hadron Collider, within the Supergravity inspired Min- 
imal Supersymmetric Standard Model. The phenomenological relevance of such re- 
actions is twofold. Firstly, they constitute a novel production mechanism of Higgs 
particles, either through a decay of a heavier (anti)squark into a lighter one or via a 
Higgs bremsstrahlung process. Secondly, their production rates are extremely sensi- 
tive to the values assumed by the five input parameters of the model, this possibly 
allowing one to put stringent constraints on the latter. After an exhaustive scan of 
the parameter space, we find that the majority of such processes could be detectable 
at high luminosity, provided tan/3 is large, tan/? ~ 30 (except in the case of iiilh and 
titg^ final states, whose detection is also possible for smaller values), that the universal 
soft Supersymmetric breaking masses are in the ranges Mq ~ 500 GeV and M1/2 ~ 220 
GeV, and that the trilinear couplings are negative, Aq < 0. We also point out some 
sizable decay signatures and discuss their Standard Model (SM) backgrounds. Finally, 
we derive compact analytical formulae of the corresponding scattering matrix elements. 



1 Introduction and motivation 



'If Supersymmetry (SUSY) exists, it will be discovered at the next generation of hadronic 
machines', has been a recurring motto so far. Indeed sooner (at the Tevatron, ^/s = 2 TeV) 
or later (at the Large Hadron Collider (LHC), y/s = 14 TeV), depending on the mass scale 
of the Higgs bosons and of the Superpartners of ordinary matter, several Supersymmetric 
'signatures' should clearly be viableQ. Typical SUSY events at hadron colliders will involve 
either the production and decay of heavy spartons, squarks and gluinos, whose foreseen mass 
range is expected to be around the TeV scale f^, or of Higgs bosons 0], primarily of the 
lightest one, for which the SUSY theory imposes a stringent mass bound of the order of the 
electroweak (EW) scale. 

However, even assuming that such a discovery will take place, there might well be little 
to learn about the fundamental dynamics of SUSY from such new events. In fact, although, 
e.g., the LHC is able to produce gluinos and squarks with masses up to 2 TeV or so and their 
detection has been shown to be feasible with rather little effort , it is much more difficult to 
determine exactly the SUSY masses involved, because in most models (i.e., those assuming 
i?-parity conservation) there are at least two missing SUSY particles in each event. Clearly, 
failing the knowledge of the SUSY mass spectrum, other typical SUSY quantities, such as 
couplings, decay rates, etc., cannot be assessed either. Needless to say, their measurement 
would be of paramount importance in order to constrain the free parameters entering the 
SUSY Lagrangian. However, by resorting to specific kinematic distributions 0, it is at least 
possible to make precision measurements of some 'combinations' of SUSY masses, but only 
in a few fortunate cases these can lead to strong constraints on the theory and its parameters. 
Besides, in minimal SUSY theories, the Higgs sector typically (i.e., at tree level) depends on 
only two such parameters, the ratio of the vacuum expectation values (VEVs) of the Higgs 
fields and one of the masses of the five physical states corresponding to the latter, as all 
others SUSY inputs enter through higher perturbative orders. Therefore, even the detection 
of a SUSY Higgs signal would carry very poor information in terms of the underlying SUSY 
model. As a matter of fact, a second question about SUSY has to legitimately be risen. 
Namely, 'Which Supersymmetric model will one discover ?' 

Thus, the key task for the Tevatron and the LHC is not only to find SUSY, but also 
to assess which model is behind it and the value of its parameters. For example, in the 
context of Supergravity (SUGRA) inspired models [0, with the minimal particle content of 
the MSSM (henceforth denoted as M-SUGRA, that we take to be the reference framework 
of our analysis) |], the dynamics of the theory can be specified by only five entries: (i) a 
universal scalar mass Mq; (ii) similarly, that for the gauginos M1/2; (hi) the universal trilinear 
breaking terms Aq (all defined at the Grand Unification scale Mgut 10 )• After the radiative 
EW symmetry breaking has taken place, two further parameters are needed to describe the 
low energy dynamics: (iv) the mentioned ratio of the VEVs of the two Higgs fields, denoted 
by tan/3 = v/v' and defined at the EW scale; and (v) a discrete parameter, sign(/i) = ±1, 
being /i the Higgsino mass term. 

Assuming universal soft breaking terms at the GUT scale, one is then able to calculate 
^For some reviews, see, e.g., jl] and [1). 
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the masses of SUSY (s)particles, their couphngs, decay rates, etc., at the EW scale, through 
the evolution of the renormalisation group equations (RGEs), the latter involving Mq, M1/2, 
Aq, tan/3 and sign(/x) as inputs. Ultimately, a comparison of such predictions with the cor- 
responding experimental measurements, as reconstructed from the actual data via dedicated 



Monte Carlo (MC) simulations [jlO[^|lT2[, should allow one to impose indirect constraints on 
the above parameters. Indeed, an additional procedure to follow in order to determine the 
latter could well be to search for the evidence of some more exotic signals of SUSY, in which, 
however, the dependence on such parameters is somewhat more manifest. 

Elementary processes of the type 



where g^'-* = t, 6, x^''^ = 1)2 and ^ = H,h, A, H^, in all possible combinations, as appropriate 
in the MSSM, serve the double purposes of: 

1. furnishing production mechanisms of Higgs bosons of the MSSM, both neutral and 
charged, in addition to the Standard Model (SM)-like channels [H; 

2. yielding production rates, for particular combinations of q^'\ x^'^ and $, strongly de- 
pendent on some of the fundamental SUSY parameters of the M-SUGRA model. 

The importance of the first point should be understood in the following terms. On the 
one hand, the detection of all neutral Higgs particles H, h and A of the theory is not certain, 
neither at the Tevatron nor even at the LHC |@]. In addition, the discovery potential of 
heavy charged Higgs bosons if ^ at both the above colliders has been proved to be extremely 
limited |jl3|. Under these circumstances, the possible existence of novel and detectable Higgs 



production channels represents a phenomenologically important result per se. (Notice that, 
for certain choices of sign(/i), Aq and tan/?, the squark-squark- Higgs coupling can become 
the largest EW coupling of the SUSY theory, even exceeding the standard Yukawa ones.) On 
the other hand, the fact that in processes (|I]) the Higgs bosons are produced in association 
with squarks via a Yukawa bremsstrahlung or in 'non-dominant' squarks decays^, implies 
that the Higgs mechanism can be probed in the sparticle sector too. In fact, other known 
production and decay mechanisms used to detect MSSM Higgs bosons mainly involve Higgs 
couplings to ordinary matter. The only exceptions are the squark loop-contributions to 
neutral Higgs boson production via gluon-gluon fusion and to Higgs boson decays through 
pairs of photons/gluons |T^, which are however swamped in both modes by the dominant 



terms involving ordinary heavy particles. 

As for the second point that we put forward, we should remind the reader of the actual 
form of the mentioned squark-squark-Higgs vertices in the MSSM (which can be found, e.g., 
in In many of these, namely when x 7^ x'; both the low-energy SUSY parameters /x 

and tan/5 enter explicitly in the Feynman rules, other than implicitly in the scalar masses. 
Furthermore, those vertices also contain A^^, the trilinear couplings at the EW scale, which 

^So that the corresponding partial widths are significantly different from the total widths, thus retaining 
the dynamics of the squark-squark-Higgs production vertices also at decay level. 
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depend critically on their common value at the GUT scale, Aq. (In the case $ = A such a 
dependence is also not affected by the mixing between the chiral, = L, R, and physical, 
X^'^ = 1,2, squark states, so that no additional SUSY mixing parameters enter the phe- 
nomenology of pseudoscalar Higgs production, this rendering the latter an ideal laboratory 
to study M-SUGRA effects [|16[). Therefore, one should expect a significant dependence of 
the production rates of the scattering processes (|I]) on tan P, Aq and /i (particularly, its sign), 
this possibly yielding a new profitable mean to constrain the underlying SUSY model. Even 
more so in the case tan /3 has previously been determined, for example, through a discovery 
in the MSSM Higgs sector. 

Concerning previous literature on the subject, we should mention that reactions of the 
type ([1|) were were first considered in Ref. Jl^ for the case gg —>■ tit\h in the so-called 
'decoupling' limit. Adopting the M-SUGRA scenario, associated production of both neutral 
and charged Higgs bosons production with squark pairs - with a special emphasis on CP- 
odd Higgs boson production - was first consider by the authors in ||16|. Furthermore, in 



Ref. |TB[ (and also jT^) light Higgs boson production in association with light top squarks 
was reanalysed in the M-SUGRA scenario at both Tevatron and LHCQ. A general consensus 



on the possible detectability at the LHC of gg titlh events emerged from Refs. 0, |19 
in the case of light top squarks and large trihnear coupling. 

We generalise here those studies, as we consider the production of all possible Higgs states, 
^ = H, h, A and H^, for a broader spectrum of their masses, in conjunction with both squark 
flavours that can have sizable couplings to Higgs particles, q = t and b, the latter taken as 
not degenerate in mass, also allowing for Higgs production in decay channels, when x x'- 
On the other hand, to simplify the simulation, we restrict ourselves to the case of gluon-gluon 
induced processes only, thus neglecting the case of quark-antiquark scatterings, which was 
instead considered in Refs. [0, |18[. This is however not restrictive. In fact, we have verified 
that at the LHC the gg contribution is around two order of magnitudes larger than the QQ 
one, in line with the findings of Refs. |]T^ |T^, well below the level of uncertainties arising in 
our computation from other sources (such as structure functions, QCD i^-factors, etc.)0. As 
for the Tevatron, we can anticipate that the production cross sections of processes (|T]) are 
generally very small (see also Ref. |18[), indeed below the level of detection over most of the 
M-SUGRA parameter space, so that we neglect further consideration of this machine here. 

It is the purpose in this paper to assess the possible relevance of points 1. and 2. in 
phenomenological studies of SUSY to be carried out at the LHC. In particular, the plan 
of the paper is as follows. In the next Section we describe how we have proceeded in our 
calculations. In Sect. |^ we illustrate the theoretical model we have resorted to in our analysis. 
Sect. ^ presents some numerical results. A brief summary and our conclusions are given in 
Sect. E. Finally, we collect in Appendix the relevant analytical formulae that we have used. 



•^The same final state but produced in e+e^ annihilations at the TeV scale, e.g., at the Next Linear 
Colhder, has been considered in Refs. 

■^As for 7, Z and s-channels, these are typically smaller by a factor of the order ©(ofem/as)^; whereas 
©(as)^ gluino interactions in t,w-channels are suppressed by the small (EW induced) mixing between light 
quarks and sbottom and stop squarks, as already remarked in ||l7|. 
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2 Calculation 



The leading-order (LO) Feynman diagrams associated to processes of the type (|l|) in the 
unitary gauge are depicted in Fig. |^. The reader can find an analytical expression of the 
corresponding matrix elements (MEs) in the Appendix. As a test of the correctness of 
our amplitudes, we have verified that they are gauge invariant by checking various Ward 
identities of the theory, both analytically and numerically. The amplitudes squared have 
then been integrated over a three-body phase space, using VEGAS and convoluted 
with gluon Parton Distribution Functions (PDFs), as provided by CTEQ(4L) [^. The latter 
constitutes our default set, taken at LO in order to be consistent with our approximation 
in calculating the scattering MEs. However, in order to estimate the systematic error due 
the gluon behaviour inside the proton, we also have resorted to other LO packages, such as 
MRS-LO(09A,10A,01A,07A) 0. Typical differences among PDFs were found to be less 



than 15-20%. The centre-of-mass (CM) energy at the partonic level, Q = V^, was the scale 
used to evaluate both the structure functions and the strong coupling constant (see next 
Section for the treatment of the latter). 

Depending on the relative value of the final state masses in (0), m^^, m-u and m$, the 

x' 

production of Higgs particles can be regarded as taking place either via a (anti)squark decay 
(if rriq^ > m-u + m$ or m-u > rriq^ + m$) or via a Higgs-strahlung (if rriq^ < m-u + m$ 

x' x' x' 

and m-'* < rriq^ + m$). In the first case, to prevent our MEs from becoming singular, 

x' 

we need to insert a finite width in the resonant (anti)squark propagators, which we have 
done by adopting the Breit-Wigner expression given in the Appendix and the appropriate 
numerical values for the widths, calculated as described in Sect. ^. Also in the second 
case, though no poles exist in the amplitudes, a finite width value has been retained in the 
propagators. Notice that we have treated the two processes on the same footing, without 
making any attempt to separate them, as for the time being we are only interested in the 
total production rates of the 2 — > 3 processes (|l|), rather than in their subsequent decay 
distributions. 



3 The theoretical model and its parameters 

In this paper we are going to display our results for squark-squark- Higgs production via 
processes (P by assuming possibly the simplest scenario in the choice of the soft SUSY 
breaking parameters at the GUT scale. That is, the so-called minimal Supergravity scenario 
or M-SUGRA inspired model, as already intimated in the Introduction. In this scenario, 
the whole dynamics of the MSSM (which contains over hundred parameters in the case of 
conserved i?-parity) at the GUT scale is reduced to the three basic inputs already introduced: 
Mq, Ml/2 and Aq. The large top Yukawa coupling then triggers the radiative EW breaking 
through the running of the soft Higgs breaking masses, from the GUT scale down to EW 
regime. From the minimisation conditions of the potential one can define the soft Higgs 
mixing parameter, and the absolute value of the Higgs mixing parameter of the SUSY 
potential, ii. The model leaves the sign(/x) and the value of tan/5 as further undetermined 
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parameters at the EW scale. 

All five M-SUGRA parameters enter into the relevant Feynman rules for the squark- 
squark-Higgs vertices, either explicitly or implicitly (through the RGEs). These can be 
written in the physical squark basis qi^2 as 

^$9292 ~ ^ q' ^^qLq'i^ ~^ ^q^q' ^'^qm'n ~ ^ q^q' ^'^qLq'^ ~ CqSqiX^q^q'^, 

-^$919^ = ~CqSqiX^qj^qi^ + SqCq/X^qj^qi^ + CqCqiX^qj^qi^ — SqSq'X^qj^qi^, 

'^*929'l = -SqCq'X^q^g'^ + CgSg' X^g^^g'^ - SgSg'X^g^g'^ + CgCg' X^g j^g'^ , (2) 

where Ql^r or ^ can in principle be any fiavour of chiral squarks. However, here we only 
focus our attention to the case of the the third generation of down and up squarks only, 
namely, sbottom and stop scalars, whose physical mass eigenstates are denoted by 61^2 and 
ti^2, respectively, the subscript 1(2) referring to the lightest (heaviest) of them. As usual, the 
Higgs fields are denoted by the generic symbol $, where ^ = H,h, A, H^. All the A<j,^ 's 
appearing in eq. (|) are function of yU,tan/5 and At^ and can be read directly from the 
Appendix of Ref. [|T^. (We ignore the case of complex /i and Aq (g = t,b) parameters by 
assuming that their phases are very small, the preferred case following the measurements of 
the Electric Dipole Moments p3|.) 

Also the left-right squark mixing angles Sg = sin6'g and Cg = cos 9q (here, q = t,b) depend 
on the M-SUGRA parameters, since they read as 

, foQ\ 2mt{At- cot (3) 

tan(26't) = , (3) 

M?^ - M?, + (i - ^)M| cos 2/3 



tan(26';, 



2'mb{Ab — yUtan/5) 



M? -M?+(-i + ^)M| cos 2/? 



with Mz the Z-boson mass and sw = sin^ 9w the sine (squared) of the Weinberg angle, 
rut and mb the top and bottom quark masses, where At and Ab are the trilinear couplings 
defined at the EW scale, while Mq^, Mqc and Mf^c are the running soft SUSY breaking 
squark masses of the third generation, for which we assume the values obtained from their 
evolution starting from a universal mass at the GUT scale equal to Mq. 

Regarding the numerical values of the M-SUGRA parameters adopted in this paper, 
we have proceeded as follows. For a start, we have set Mq = M1/2 = 150 GeV. Such 
rather low values for the universal masses come as natural first choice, if one is interested 
in detecting processes of the type (|I|). For two simple reasons. On the one hand, these two 
quantities determine the actual rrig^, nig' ^ and m<i, values entering processes (|l]), through their 
intervention in the RGEs, in such a way that small values of Mq and M1/2 at the GUT scale 
convert into a rather light squark and Higgs mass spectrum at the EW scale. On the other 
hand, being 2 — > 3 body processes, a strong suppression from the phase space would arise 
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in squark-squark-Higgs production if the masses in the final state were too large^. For the 
above choice, the M-SUGRA model predicts squark and heavy Higgs masses in the region 
of 80-450 GeV (as we shall see in more detail below), so that the latter can in principle 
materialise at LHC energies (further recall that the light Higgs mass is bound to be below 
130 GeV). 

Then, we have varied the trilinear soft Supersymmetry breaking parameter Aq in a region 
where it changes its sign, e.g., (—300, +300) GeV, while we spanned the tan/? value between 
2 and 40. As for /i, whereas in our model its magnitude is constrained, its sign is not. Thus, 
in all generality, we have explored both the possibilities sign(yu) = ±1. 

As a further step of our analysis, we have then come back to Mq and M1/2 and change 
them, while maintaining tan (3, Aq and sign(yu) fixed at some specific values. We have done 
so only for those processes that we had already identified to have not only a large cross 
section, but also a strong dependence on one or more of these three M-SUGRA parameters. 

Given the strong phase space suppression induced by the consequent increase of mg^, 
rriq' and m<j, in the final states, we will cautiously maintain the universal scalar and gaugino 

x' 

masses below 250-300 GeV (at least at first). However, the reader should not assume that 
this is a necessary condition to the experimental detection of processes of the type (|l]). In 
fact, this need not be true, as we shall show that even for Mq values as large as 500 GeV 
one can find sizable squark-squark-Higgs production cross sections for M1/2 up to 200-250 
GeV, as long as Ao is strongly negative, tan/? ~ 30 and sign(/i) < 0. 

Such unexpected behaviours are strongly driven by the intervention in the production 
rates of tan/?, Aq and sign(/i), through the trilinear scalar couplings \<s>q^q' niore than by 
the actual values of rrig^, rriq' ^ and m$. This is evident by a mere look at the standard 



Feynman rules, as can be found, e.g., in Ref. [|T5[. We will make our concern in this paper 



that of guiding the reader through such delicate interplay between masses and couplings, 
by explicitely writing down the expression of the relevant vertices in those parameter space 
domains where such complicated phenomenology manifests itself. 

Starting from the five M-SUGRA parameters Mq, M1/2, Aq, tan/3 and sign(/i), we have 
generated the spectrum of masses, widths, couplings and mixings relative to squarks and 
Higgs particles entering reactions (|l|) by running the ISASUGRA/ISASUSY programs contained 
in the latest release of the package ISA JET ||10|, version 7.40. The default value of the top 
mass we used was 175 GeV. Note that also typical EW parameters, such as and sin^ 9w, 
were taken from this program, as they enter the RGBs of the SUSY theory. Concerning 
the value of the strong coupling constant, as, entering the production processes (p, we have 
proceeded as follows. By using as inputs the extracted value of ag at the Mz scale, we evolve 
it up to any scale Q by making use of the two-loop renormalisation group equations and by 
taking into account all the low-energy threshold effects from the various SUSY masses by 



means of the theta function approximation, as discussed in Ref. [26 



^The additional depletion coming from the gluon PDFs, which would be probed at much higher values of 
(of the order of the rest masses or more), where they are naturally smaller, would in part be compensated 

by the rise of the quark-antiquark initiated subprocesses: i.e., Q + Q ^ qx+ Q^' + ^, where, again, one has 

that qC) = t, b, x^'^ = 1, 2 and $ = H, h, A, H^. 
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Finally, notice that in scanning over the M-SUGRA parameter space, one should make 
sure that the values generated for the Higgs boson and sparticle masses are in accordance 
with current experimental bounds. Signs of the sort "x" or shaded areas appearing in our 
figures in forthcoming Sect. ^ will correspond to already prohibited areas for the parameter 
space of our SUSY model. We nonetheless leave them for illustrative purposes, in order to 
visualise the typical impact of present and future experimental bounds on the phenomenology 
of our reactions. For example, the M-SUGRA points individuated by the combinations 
Mq = Ml/2 = 130 — 150 GeV, tan/3 = 2, Aq = GeV and sign(/i) = — , used in some of 
the tables and figures in the next Section, contradict the limits on the lightest Higgs boson 



mass from direct searches as they yield nih = 72 — 80 GeV. We will discuss the 

experimental bounds in more detail in the next Section. 

As for theoretical constrains, these arise from two sources: namely, the absence of charge 
and colour breaking minima and that of large contributions to the EW observables that are 
measured with high precision at LEP. The former is avoided when the following inequalities 
hold tree level E^: 





< 


3 


Al 


< 
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Al 


< 
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2 , 2 I 2 I 2 



is 



Al < 3(m|^+m|. + /i2 + m^J, (4) 

where all masses appearing in eq. (^) are the soft Supersymmetry breaking masses except 
the Higgsino mixing parameter ii. When Aq is below 1 TeV, as is the case in our analysis, 
the above constraints are always satisfied even for very light squark masses. As for the 
contributions to the EW observables, we have found the region 150 GeV ~ Mo,Mi/2 ~ 500 
GeV, 2 ~ tan/3 ~ 40 and |Ao| ~ 900 GeV covered by our analysis in agreement with the most 
recent measurements of the 'effective' weak mixing angle, sin^ 9^,^ = 0.2321 ± 0.010, and of 
the mass, Mw = 80.388 ± 0.063 GeV In particular, notice that the M-SUGRA 
prediction for sin^ 9es decreases for a lighter mass spectrum while it becomes constant in the 
heavy mass region 



4 Results 

We begin this Section by analysing all reactions (||) in the low mass regime, i.e., that induced 
by values of Mq and M1/2 below 250-300 GeV. This is done in Subsect. |]T. In this scenario. 



we will first present and discuss, for future reference, the values of squark and Higgs masses 
resulting from the RGE evolution: in [4.1.1| and ^4.1.2| , respectively. Then, we will move on 



to considering the production of, in turn: neutral CP-even (in [4.1.3| ), CP-odd (in [4.1. 4|) and 
charged (in [4.1.5| ) Higgs bosons. Possible decay signatures of the latter will be analysed in 
4.1.6| . Finally, Subsect. will pin-point those unusual cases discussed above, in which the 



suppression from very heavy scalar masses in the final states of reactions (|l|) can be overcome 
by strong vertex effects, yielding in the end sizable cross sections. 
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4.1 Light mass spectrum 



Once fixed Mq = M1/2 = 150 GeV, one obtains the (sbottom and stop) squark and Higgs 
masses reported in Figs. ^ and ^ respectively, depending on 2 ~ tan /3 ~ 35 and for Aq = 
—300, 0, +300 GeV. Tlie two possible options for the sign of /x are also contemplated. 

Far from willing to discuss exhaustively the dependence of rriq^ , m^' ^ and m$ upon the 
five M-SUGRA parameters, we limit ourselves here to spotting in Figs. ^|-^ some interesting 
trends, that will affect the overall behaviour of the squark-squark-Higgs cross sections that 
we will be treating below. For a more complete overview, see, e.g., Refs. [B^ . 



4.1.1 Squark masses 

The four squark flavours of the first and second generation [q = u, d, c, s) with left- and right- 
handed components are all nearly degenerate in mass and the latter is given approximately 
by the following formula: 



m. 



Ml^^Ml {q = u,d,~c,~s). (5) 



g — Y ""0 I ""•'-'1/2 

For our choice of Mq and M1/2, one gets rriq ^ 400 GeV. The light squark flavours are not 
our concern in processes (0), though they may well enter some of the decay chains of the 
other (pseudo) scalar particles produced. 

The two squark flavours of the third generation must be treated differently because the 
off-diagonal entries of their mass matrices can be large, owing to the strength of the Yukawa 
couplings of the corresponding quarks (in this respect, notice that the bottom one becomes 
comparable to that of the top in the large tan/3 region |T^). It thus follows that the mass 
eigenstates ti, t2, bi and 62 are all different and generally smaller than rriq. Among these, ti 
is most often significantly lighter than all the other stop and sbottom states. At large tan j3, 
the same happens to the bi mass eigenstate, as large values of tan/3 correspond to smaller 
sbottom masses, so that one eventually gets that m^^ ~ m^^ . 

Variation of the trilinear couplings can also cause significant differences between the light 
stop and sbottom masses. Finally, the sign of the Higgsino mass term plays an important 
role when tan/5 is either small or large, for the cases w^t^2 ^^"^ '^fci2' respectively. 

Experimental limits on the squark masses come from searches at Tevatron and LEP2. 
The most stringent bound on the ti mass comes from the hadron collider [35|: in absence 
of mixing, values of m^^ < 120 GeV are excluded foi|^ ^lsp = ^ '^^ GeV. D0 exclude 
values of m^^ below 85 GeV for m^o < 47 GeV |Q and ALEPH do over the region m^^ < 83 



GeV for any value of LSP mass In addition, CDF exclude masses for the lightest 



top squark up to 120 GeV when the LSP is m^o < 50 GeV Finally, D0, using data 
corresponding to an integrated luminosity of 79 pb^^, contradicts all models with 'm'-_^^_^ < 
250 GeV for tan/? ~ 2, = GeV and fi < (in scenarios with equal squark and 
gluino masses the limit goes up to fn-_^^_^ < 260 GeV). 

^Recall that in M-SUGRA the Lightest Supersymmetric Particle (LSP) is the lightest neutralino, Xi- 
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4.1.2 Higgs boson masses 



The mass of the hghtest Higgs boson is here constrained to be less than 115 GeV and it 
appears to exhibit constant values over the region 8 ~ tan/5 ~ 35, for a given combination of 
Ao and sign(/i). To change either the trilinear coupling or the sign of the Higgsino mass has 
the net effect of scaling rrih, lower with increasing Aq and higher for positive fi. 

As for the masses of the other Higgs bosons, they are nearly degenerate. They vary 
between 180 GeV (when tan/5 is large) and 400 GeV (when tan/3 is small). The dependence 
on sign(/i) is generally negligible, whereas the one on Aq is very strong in the intermediate 
tan f3 regime. One thing worth noticing here is the existence of a point where all three Higgs 
masses mH,rnA and converge, regardless of the values of Aq and sign(/i). This occurs 
for tan/3 f» 30, where 

rriH = rriA = fnfj± ^ 200 GeV. (6) 

As for experimental bounds, LEP experiments have combined their results from data 
taken at CM energies from 91 to 183 GeV to place lower bounds on the masses of the 
light (m/i) and pseudoscalar (m^) Higgs bosons, of 78.8 and 79.1 GeV, respectively [^. In 
addition, they exclude the range 0.8 < tan (3 < 2.1 for minimal stop mixing and rrit = 175 
GeV. Also, D0 have recently removed at 95% confidence level the intervals tan/3 < 0.97 



and tan^ > 40.9 for Mh± = 60 GeV and a{tt) = 5.5 pb (again, with mt = 175 GeV). 
However, the limits become less stringent with increasing Mij±\ e.g., for Mh± > 124 GeV 
(as is the case here) the available angular range is 0.3 < tan/3 < 150. 

4.1.3 CP-even Higgs boson production 



We present in Figs. ^-|^ the tan /3 dependence of the production cross sections of the two 
neutral scalar Higgs bosons, h and H, respectively, in association with any possible com- 
bination of squarks of the third generation. Again, we parametrise the dependence upon 
Aq by adopting for the latter the discrete values of —300, and +300 GeV and we choose 
sign(/i) = ±1. 

Assuming an integrated luminosity of 100 inverse femtobarns over a twelve month period 
of running at the LHC (i.e., J Cdt = 100 fb~^), one can realistically hope for the detection 
of squark-squark-Higgs processes only if the production cross section is above 1 fb or so. In 



fact, we shall see in Sect. |4.1.6| how typical decay fractions of clean signatures range at the 
level of 10% or below (see Subsect. [1.1. 6| later on). 

Under this assumption, one immediately sees that there are several production channels 
of CP-even Higgs particles which could be observed, over a large part of the M-SUGRA 
parameter space considered here. Primarily, those involving the lightest stop squark, ti, 
particularly if also the lightest Higgs state is involved, but not only. 

For the case of h production, there exists an approximate hierarchy of cross sections 
which can possibly be detected: 

a{tit\h) ^ a{titlh) > a(6i6*/i). (7) 

The cases 62^2^, ^2^2^ and hib^h never have significantly large rates. 
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In the case of titlh states, those with largest production rates, one obtains 

a{gg htit\) > 20 fb, (8) 

for every combination of tan/5, Aq and sign(/i), in the case of Mq = M1/2 = 150 GeV, thus 
a sort a lower limit over a representative portion of the low mass regime of the M-SUGRA 
scenario. Moreover, the largest production rate for this final state (compatible with the 
current experimental constraints) is obtained in the small tan/5 ~ 2 region and for the 
combinations Aq ~ —300 and /i = +, for which 

a{titlh) > 200 fb, (9) 

corresponding to more than 20000 events per year running of the LHC 

The dominance of the production channel involving both the lightest squarks and Higgs 
boson, above all other mechanisms (|T]), was foreseeable. The reason is rather simple, in fact, 
twofold. On the one hand, the sum of the rest masses in the final states yields the smallest 
possible values, thus enhancing the volume of the three-body phase space, relatively to 
any other squark-squark-Higgs combination. On the other hand, the cross section is also 
significantly enhanced when the trilinear coupling assumes negative values. From the 

analysis of the RGBs we find that At 400(-250) GeV when Aq = -300(300) GeV. Now, 

the coupling of the light Higgs boson to top squarks is driven by At when the latter takes 
on large values. More specifically, the corresponding vertex reads as (recall that we are in 
the kinematic limit mx ~ itlh ~ ^ ^ Mz'- see Fig. |]) 



htiti 



igwMz 

Cw 



f r 1 2 

< — cos 2/3 cos^ 9t H — cos 2/5 cos 26*^ 
L . 2 3 



mf nit sin 29t 



M| 2M| 



(A + /icot/5)} , (10) 



(here and in the following, g'^ = Airaej^/s'^, cw = ^1 — and = M|(l — s^)) where 
for large At (note that sin 9t is maximal in such a case) the coupling goes like XhiJi ^ ^m^- 
As a consequence, because of the presence of the trilinear term At, the coupling of the light 
Higgs boson to light stop squarks could be much larger compared to that to the top quark, 
which behaves like Xhtt oc jj^, as already recognised in [^, 0. Over the parameter space 
that we have chosen here, the cross section of gg t\t\h can be either larger than or of the 
same order as that of gg tih |TB|. 



Thus, the subprocess gg tit\h can well boast the status of an additional discovery 
mechanism of the lightest scalar Higgs boson of the MSSM, as remarked in Refs. ||l^ |18[| . 
(This is true also in non M-SUGRA models, where however one still has that mn ~ ~ 



mH± 3> nih IT^, |T8l-) In this respect, the reader should further notice the stability of its 



production cross section against variations of tan/5 (see also, e.g.. Fig. 4 in ||TH|)[]. This 

^The tan/3 dependence of a{t\t\h) at low values of such a parameter is mainly a phase space effect, as 
it can be deduced by comparing Fig. || to Fig. ^. In addition, in the low tan/3 domain, there are residual 
effects onto the production rates induced by the term /i cot (3 arising from the off-diagonal elements of the 
squark mass matrices and affecting the A^^^^^ vertex. 
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proves to be a crucial point, as it is possible that one will have no narrow hints about the 
actual value of this crucial parameter of the Higgs sector even after Run 2 at Tevatron 
(unless, of course, the lightest Higgs boson is discovered there !) In other terms, iiilh 
would always be present at fixed rate at the LHC, no matter whether tan (3 is large or small. 
Similar arguments can be put forward concerning the sign(/i) dependence. 

Some care must instead be exercised with respect to the Aq dependence. In fact, to vary 
the universal trilinear coupling between, e.g., —300 and +300 depletes the cross section by 
a factor of about seven, as shown in Fig. ^. For even larger differences, say, between —500 
and +500 (not shown here for reasons of space), the ratio between the cross sections become 
as big as 30 ! Not surprising ly then, Ref. [||] focused on the choice Aq = -2000 GeV0 
(and sign(yu) = +). Far from regarding this dependence as a shortcoming of titlh events in 
helping in the quest for the so far elusive lightest Higgs boson of the MSSM (in fact, even for 
very large and positive Aq values we found the cross section well above 1 fb), this example 
allows us to enlighten that other aspect of squark-squark-Higgs production that we have 
mentioned in the Introduction: i.e., its potential in pinning down some of the fundamental 
parameters of the M-SUGRA model. Needless to say, variations of the cross section with 
Aq as large as those mentioned above are well beyond the various sources of uncertainties 
on the production rates (other than the theoretical ones related to the PDFs and the effect 
of higher-order QCD corrections also the experimental ones in their determination). To 
measure a production rate of titlh events much larger than about 50 fb (the value obtained 
in correspondence of Aq = 0), in some specific decay channel, would unambiguously imply 
that the universal trilinear coupling at the GUT scale is negative. 

As already mentioned, very little could be learn about the actual values of tan/? from 
this specific process. However, if the latter is known beforehand to be around 2 or so, one 
could use this information to constrain sign(/i). In fact, for Aq = —300 GeV, one would get 
that 

aiiJlh) ^ 200 fb ^ /i = + (11) 
a{iiilh) ^ 50 fb ^ ^ = - . (12) 



Let us proceed in this spirit to see whether other channels can be of some help in con- 
straining the M-SUGRA model. Following the list of detectable h production cross section 
given in (^, we find the ^1^2^ final state |jT6|. This is not surprising either. In fact, the 



relevant coupling behaves like (again, X = H, A, if^) 



Xr^x»'^'z ^ 1 COS 2/? - - -s'^ sin 



(A + /icot/?) , (13) 



rrit cos 29i 
2M| 

becoming very large when cos2^j and {At + yucot/?) reach their allowed maximum values. 

^Note that in this ^0 region the M-SUGRA scenario clashes against the constraints from the charge and 
colour breaking minima, i.e., eq. (^. This is the reason why we prefer to display our results in a rather more 
conservative range, i.e., \Aq \ ~ 1 TeV. 
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From Fig. one can see that this happens in the region of small tan /5, negative sign of /i 
and Aq < -300 GeV. 

The intriguing aspect here, which was largely missing in the case in which both squarks 
were the lightest, is that one could impose severe constraints on the sign of the Higgsino 
mass term, other than on Aq. In fact, in the detectable region, the curves corresponding 
to sign(/i) = — (higher) and sign(/i) = + (lower) depart considerably. For example, for 
relatively small tan/? values, say 4, the ratios as obtained by dividing the cross sections 
corresponding to negative /i's by those for positive Higgsino masses are quite large indeed: 
about 7(5) [2] when Aq = —300(0) [+300] GeV. At even lower tan (3, say, equal to 2, one 
symbolically has: 

a{iii;h) ^ 300 fb ^ = - (14) 
aiiJlh) ^ 2ih ^ fi = + , (15) 

(e.g., for Aq = —300 GeV). Luckily enough here, where the solid and dot-dashed curves 
start getting closer (for tan/3 ~ 15 — 20) is precisely when the cross section is no longer 
observable. However, what just said makes the point that tan (3 ought to be known rather 
accurately from some previous measurements, if one wants to constrain the other M-SUGRA 
parameters by studying the production of the lightest Higgs scalar of the theory produced 
in association with both stop mass eigenstates. 

In our list of observable h cross sections biblh comes next. Here the potential is somewhat 
complementary to the two above cases, in the sense that not to find any pairs of sbottom 
squarks of the type bi produced in association with an h scalar once tan (3 is already known 
to be large could have powerful consequences on the viability of M-SUGRA as the underlying 
model of SUSY. To be specific, notice how the six curves corresponding to all the possible 
combinations of the parameters Aq = —300, 0, +300 GeV and sign(/i) = ± lie within a factor 
from 2 to 4 in cross section, in correspondence of tan/5 = 20 and 35, respectively. Even the 
cases of Aq = ±500 GeV do not depart significantly from the central curve for = 0, in the 
above tan/? region. Unfortunately, contrary to the case of titg^ production, here the most 
interesting region is presumably below detection level. In fact, for tan/? quite low, a huge 
portion of M-SUGRA parameter plane collapses into a narrow stripe, as the various curves 
tend to overlap, all being contained within a factor as small as 1.5 (e.g., at tan/3 = 2, also 
including the cases |v4o| = 500 GeV, not shown in the figure). 

As for H production, one identifies as possible candidates for detection the following 
cases: 

a(biblH) ~ a{iiilH) > aihblH) ~ a{iii;H). (16) 

The remaining two combinations, i.e., bib^H and t2^2-^; yield cross sections that are hope- 
lessly small. 

Also some of the detectable H production processes can have a significant impact in 
aiding the determination of the M-SUGRA parameters, most notably those yielding the 
final states bib^H and 62&2-f^- Here, if tan/? is known to be, say, 35 or so, a detection of 
either the former or the latter by the thousand or hundred, respectively, would imply that 
Aq is most certainly negative, since production rates corresponding to Aq values larger than 
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zero are about a factor of 5 and 3 smaller (rather irrespectively of sign(/i)). Somewhat less 
discrimination power between positive and negative Aq values have tit\H and tit^H events, 
over the same (large) tan/3 region as above. The most interesting case would have been 
^2^2-^^, as the collapse of the M-SUGRA model that we already noticed in the case of 
final states is here even more striking, over a more significant tan j3 region. Unluckily enough, 
the corresponding production cross section never exceeds the femtobarn level. 

A general remark in now in order, concerning the strength of the coupling of sbottom 
squarks to neutral CP-even Higgs bosons. The monotonic growth of the production rates of 
sbottom squark processes with increasing tan /?, as opposed to a much milder dependence 
of the stop ones, has a simple explanation. For example, the coupling A^j^j^ is driven by 
the term ^^Bdi^ and its size becomes large for large tan/5. Another reason for an extra 
enhancement of the sbottom production rates comes from the phase space available to the 
final states, as both m^^ and m^^ decrease very fast when tan (3 gets large, whereas this is 
much less the case for m^^ and mi^ (see Fig. There is however a point, for the case of 
reactions involving the two sbottom mass eigenstates at once (i.e., 6162^ and hib^^H), in which 
the production cross sections vanish altogether, somewhere in the vicinity of tan /5 = 34 — 36 
(the zero for bib^H is actually beyond the tan/3 range plotted), the exact value depending 
upon Aq and sign(yu). This is clearly induced by the Xf^i^i^ and ^Hbib2 vertices and their 
typical oc (/i — Afetan/?) behaviour, when -C |Afetan/3| and Ab changes its sign. 

Before closing this Section, we investigate the residual dependence of CP-even Higgs 
boson production in association with sbottoms and stops on the input values of Mq and 
Ml/2, when they are allowed to deviate from their common default of 150 GeV assumed so 
far. For illustrative purposes, we do so by adopting two discrete values of tan/3, 2 and 40, 
and choosing the combination Aq = and negative /i. Anyhow, though not shown, we have 
verified that a similar pattern to the one that we will outline below can be recognised also 
for the case of finite (positive and negative) values of Aq and positive Higgsino masses as 
well. For this exercise, we focus our attention only to the dominant production cross sections 
in either case, that is, tit\h, tit2h and hib\h for light (see Tab. 0) and tit\H and hib\H for 
heavy (see Tab. ^ Higgs bosons. 

We obtain that most of the cross sections with the light Higgs particle involved decrease 
when either or both the parameters Mq and M1/2 increase. In this respect, however, it is 
well worth noticing that the total cross section for tii^h production acquires a large statistic 
significance in the higher mass regime and maintains it even at the very upper end of it 
(some 5000 events/year can be produced at the LHC via this mode if tan/? = 2, Mq = 300 
GeV and M1/2 = 250 GeV). In this area of the M-SUGRA parameter space, tif^h events 
are much more numerous than tii\h ones, the other way round with respect to the low mass 
combination Mq = M1/2 = 130 GeV, despite of the more massive final state. (However, this 
is only true al low tan /3.) In fact, for Mq = 300 GeV and M1/2 = 250 GeV, the squark masses 
are mi^ = 472 GeV, m^^ =591 GeV, m^^ = 569 GeV and 171^^ = 623 GeV. The inversion 
of hierarchy between the two cross sections is induced by the onset of the ^2 — ^ iih decay 
channel at large Mq and M1/2 values, as m^^ ~ m^^^ + ruh, whose resonance enhancement 
in the 2 — 2 process gg ^2^2 overcomes both the inner phase space depletion and the 
strength of the Higgs-strahlung emission in gg — > iiil events. 
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Table 1: The variation of the most significant cross sections of processes gg Qx^x'^ with 
Mq, Ml/2 and tan/3. The other M-SUGRA parameters are fixed as follows: Aq = GeV and 
sign(/i) = -. 

As for the case biblh (and similarly for 62^2^; not shown in the table), there is no inversion 
of tendency here, in the interplay with the light stop channel, as production rates are strongly 
dominated by the fact that m^^ , m^^ ^ m^-^ . They both are very much suppressed. On similar 
grounds, one can argue about the smallness of t2t2h. Finally, being m^_^ ~ m^^ + in most 

part of the (Mo,Mi/2) plane that we have spanned, the 62^1^ final state never stands up 
either as quantitatively interesting. 

The results for the mass dependence of the production rates for the heavy (CP-even) 
Higgs boson have a simpler pattern. The all of their production phenomenology is governed 
by the fact that over the (Mo,Mi/2) regions considered here one never finds the kinematic 
configuration m^^^^^ ~ ''^t{b) +'^-H'- That is, no production and decay channel can onset, and 
the hierarchy of cross sections already seen for Mq = M1/2 = 150 GeV replicates unaltered 
in most cases, mainly governed by the size of the rest masses in the final state. Here cross 
sections remain sizable only if neither Mq nor M1/2 exceed the value 150-200 GeV (for large 
tan P, of course, see Fig. ^. Even in such cases, though, presumably no more than a handful 
of events can be selected in most decay channels. 



4.1.4 CP-odd Higgs boson production 



As discussed to some lenght in Ref. |TB|, pseudoscalar Higgs boson production in association 
with sbottom and stop squarks of different mass (the only possible combination in absence 
of CP- violating phases in fi and Ag, with q = t,b), can boast a special attractiveness because 
of the absence of mixing terms in the relevant squark-squark-Higgs couplings. By making 
use of eq. (H) and recalling that if one reverts the chirality fiow in the vertex XAgLQR the 



corresponding Feynman rule changes its sign |15], one finds that those vertices reduce to 



9wm 
"2M, 
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(H- Ab tan P) . (17) 
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Table 2: The variation of the most significant cross sections (in pb) of processes gg q^q^,H 
with Mq, Ml/2 and tan/3. The other M-SUGRA parameters are fixed as follows: Aq = 
GeV and sign(/i) = — . 

These are precisely the couplings entering the two processes of the type (0) inducing the 
final states gi(?2^, where q = t,b. 

From this point of view, it is then clear the potential of squark and pseudoscalar Higgs 
production in constraining the input values of all five M-SUGRA parameters. In other 
terms, to trace back (more technically, to fit) the shape of the cross sections (if not of some 
differential distributions) in terms of the tan (3, Ag (with q = t,b) and /i parameters entering 
eq. {\rf\j is presumably a much simpler job then doing the same by using the more involved 
expressions in eq. (^, unless one exploits some asymptotic regime in either tan/3, Ag (with 
q = t,b) and/or fi in which the latter reduce to the former. It is under this perspective that 



we looked at the case of A production in our Ref. ||16 |. 

Rather than repeating the all discussion carried out there, we summarise here the salient 
findings of |T^, referring the reader to that paper for specific details. The production cross 



sections can be found in Fig. ^. For tan P below 20 or so, the rates for pseudoscalar Higgs 
boson production are presumably too poor to be of great experimental help. Furthermore, in 
the high tan (3 regime, pseudoscalar Higgs boson production is in general less effective than 
other channels in constraining the sign of the Higgsino mass term: compare the overlapping 
for the solid and dot-dashed curves (for each Aq) in the detectable regions of ti^A and 6162^ 
production to the splitting occurring in, e.g., tit\H. This, as far as it concerns the flaws. 

As for the advantages, we would like to draw the attention of the reader to the fact that 
reactions (|l|) with CP-odd Higgs bosons in the final state are quite sensitive to tan /3. The 
simple form of the expressions for \Aq1q2 il = ^) eq. (|1^ allows one to straightforwardly 
interpret the variation of the pseudoscalar rates with this parameter, namely, the steep rise 
at high values of the latter. This can in fact be understood as follows. For large tan j3, the 
vertex couplings of eq. (|T7|) can be rewritten in the approximate form 



that is, the coupling which is associated with the sbottom pair is proportional to tan (3, so 
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that, eventually, the total 6162^ cross section will grow with tan^ j3 while the coupling related 
to the stop pair will assume constant values. In the latter, the enhancement of the ^1^2^ 
cross section with tan jS is rather a phase space effect since, as tan (3 increases, the CP-odd 
Higgs boson mass decreases considerably (the squark masses changing much less instead), 
as we can see from Fig. ^. Of course, the same remains valid in the former case as well, so 
that our figure indicates a clear order in the size of the cross sections, a{biblA) ~ aitit^A), 
at large tan f3. 

But, let us now turn our attention to another peculiar dependence of the production 
rates of tit^^A and 6162^: the one on the common trilinear coupling Aq. Pretty much along 
the same lines as for the combinations tit\h, tit2h and hih\H one can make the case that the 
sensitivity to Aq of the A production cross sections offers the chance of constraining, possibly 
the sign, and hopefully the magnitude, of this fundamental M-SUGRA parameter. This is 
presumably the best attribute of tif^A and 6162^, under the assumption already made in 
few instances that the determination of tan (3 could come first from studies in the pure Higgs 
sector. Putting down some numbers in this respect, one may invoke the following scenario, 
if tan (3 is, say, larger than 32: 

a{titlA) > 10 fb ^ < -300 (19) 

a{hihlA) I 2 fb ^ Ao < -300, (20) 

quite independently of sign(/i). Conversely, the non-observation of pseudoscalar Higgs events 
in those regimes would imply most likely a positive Aq value. 

As for peculiar trends in Fig. ^ it is worth mentioning (though we have not shown it 
here, as the reader can refer to |T^) that the cross section for sbottom production vanishes 
too, at some (large) value of tan (3, as it did for the case of CP-even Higgs production. As 



the reader can appreciate in Fig. 2 of Ref. [16|, contrary to the case of 6162^ and bib2H final 



states, this however happens in 6162^ events only for positive and large values of Aq (and 
both sign(/i) = ±). This is another consequence of the different nature of the couplings (|^) 
to squarks of CP-even versus CP-odd neutral Higgs bosons. Though we have failed to find 
a point where this disappearance of 6162^ events for positive Aq values corresponds to the 
survival of a detectable cross section for negative Aq^s (at fixed tan/?), such matter would 
presumably deserve further investigation in the future. 

Before closing this Section, we study the dependence of pseudoscalar Higgs boson produc- 
tion in association with stop and sbottom squarks on the last two M-SUGRA independent 
parameters, Mq and M1/2: see Tab. ^. The main effect of changing the latter is onto the 
masses of the final state scalars, through the phase space volume as well as via propagator 
effects in the scattering amplitudes. (In fact, no decay channel of the heavier stop or sbottom 
into the lighter one ever opens, at least for the values of Mq and M1/2 that we had looked 
at.) In other terms, to increase one or the other depletes both (7(^1^2^) and (7(6162^) quite 
strongly, simply because the values of all rrig^'s, rriq ,'s and m<i,'s get larger. For example, 
assuming tan (3 = 40: at Mq = M1/2 = 130 GeV, one has m^-^ = 248 GeV, m^-^ = 388 GeV, 
= 256 GeV, m^^ = 340 GeV and tua = 120; whereas at Mq = 300 GeV and M1/2 = 250 
GeV, the figures are m^^ = 461 GeV, m^^ = 611 GeV, m^^ = 510 GeV, m^^ = 591 GeV and 
uiA = 292 GeV. In practice, the table shows that only rather light Mq and M1/2 masses (say. 
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Table 3: The variation of the most significant cross sections of processes gg — > q^^q^A with 
Mq, Ml/2 and tan/3. The other M-SUGRA parameters are fixed as follows: Aq = GeV and 
sign(/i) = -. 

below 200 and 150 GeV, respectively) would possibly allow for pseudoscalar production to 
be detectable at the LHC, and only at large tan/3 [0. 

4.1.5 Charged Higgs bosons production 

To have an additional source of charged Higgs bosons at the LHC, especially with masses 
larger than the top mass rrit, would be very helpful from an experimental point of view. In 
fact, it is well known the difficulty of detecting charged Higgs scalars in that mass regime, not 
only because of a dominant decay signature which suffers from very large QCD background 
(i.e., ^ tb bbW~^ bbjj), but also because the production mechanisms are not 
many and with not very large rates [jl3|. Unfortunately, as it turns out from Fig. |^, typical 
production cross sections of scalars in association with sbottom and stop pairs rarely 
exceed 10 fb. These rates compare rather poorly with other mechanisms |T^, for the same 
choice of . Thus, there is little to gain in exploiting processes (|l|) as discovery channels 
of charged Higgs bosons. 

Furthermore, their dependence on tan /3, Aq and sign(/i) replicates many of the tendencies 
already individuated in neutral Higgs channels, for which the production cross sections are 
much larger. Adding the fact that typical decay channels of the latter (e.g., in photon pairs) 
are much cleaner in the hadronic environment of the LHC than those of the former, one 
would quite rightly conclude that the potential oiiiblH~, tifog-f^", bii^^H^ and 62^2-^^ final 
states in constraining the M-SUGRA parameter space is rather poor. 

Nonetheless, it is worth recognising some of the typical trends of the production cross 
sections, for the sake of future reference. Let alone the last two combinations, for which the 
final state masses are too heavy to be produced at detectable rate, we have a quick look at 
the first two cases, which can indeed have cross sections significantly above 1 fb, at least in 
the large tan/? region. This enhancement has a twofold explanation. Firstly, phase space 
effects, as for large tan/5 all scalar masses (except m^-^ and rrih) get smaller: see Figs. ^|-^ 
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Table 4: The variation of the most significant cross sections of processes gg — > q^q^,H^ with 
Mo, Ml/2 and tan/5. The other M-SUGRA parameters are fixed as follows: Aq = GeV and 
sign(/i) = -. 



Secondly, terms in their couplings proportional to mi,Ai,tan(3 are dominant for most of the 
possible Aq and sign(/i) combinations. Finally, notice also in the case of the tih\H^ final 
state the vanishing of the cross section, this time at somewhat lower values of tan (3 than in 
the case of the neutral Higgs bosons. 

As for the Mq and M1/2 dependence, this is again realised through the phase space and 
the propagators, as there is no significant enhancement from resonant decays. In practice, 
only if tan (3 is extremely large and both the universal scalar and gaugino masses are below 
200 GeV, the two cross sections for tih\H~ and tib^H' remain detectable (indeed, those 
containing the hghtest stop squark): see Tab. ^ 



4.1.6 Decay signatures 

So far we have only discussed production cross sections for processes of the form (Q) and made 
no considerations about possible decay channels and relative branching fractions of either 
squarks or Higgs bosons. Another related aspect is the typical kinematics of the signals, as it 
appears in the detectors, and the size of the possible backgrounds. Furthermore, the reader 
should appreciate how all channels entering processes (|I]) are intertwined, in the sense that 
any of these can act as a background to all others. 

It is the aim of this Section that of indicating some possible decay signatures of the 
most relevant squark- squark- Higgs processes, in which they show both large rates and their 
kinematics is such that they can hopefully be disentangled at the LHC. In doing so, we 
distinguish between a small (see 4.1.6| .l) and large (see 4.1.6| .2) tan/3 regime, as we have 



shown that such a parameter is crucial in determining the actual size of the production 
rates. As representative choice of the universal masses we adopt the combination with 
lowest values among those discussed in the previous Subsections, i.e., Mq = M1/2 = 130 
GeV, further setting Aq=0 and sign(/i) negative. 
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4.1.6 



1 Small tan/3 regime 



For small tan /?'s the only relevant processes are titlh and titg/^ production. A possible decay 
signature for the first case is the one contemplated in Refs. [l^, [l^. That is, ti xt^ — > 
W^b plus missing energy for the light stop and /i — > 77 for the light Higgs boson, with the 
decaying leptonically and/or hadronically. The final topology would be the same as in 
tih, with the only difference that for stop squark events there is a large amount of missing 
energy. 

Since another option to tag the lightest MSSM Higgs boson at the LHC is to use the 
more messy but dominant decay channel into bb pairs (as opposed to exploiting the cleaner 
but suppressed 77 mode) |^ another possible decay sequence could be the following: 

h il h 

\. J, 

Xt + b Xl +b b + b 

I I 

q + q' + X^ + + X°i 

in which qq' = ud, cs and £ = e, /i. Considering also the charge conjugated XiXi decays, the 
final signature would then be '2 jets + 46 +i^ + Ej^iss\ where the missing energy /momentum 
is not only due to the two Xi's but also to the neutrinos. 

The total branching ratio (BR) of such a decay sequence is, for tan /3 = 3 and according 
to Tab. ^ approximately 2.5%. The production cross section at the same tan/3 value is 
about 72 fb, so that about 176 events per year would survive. One may further assume 
a reduction factor of about 0.25 because of the overall efficiency to tag four 6-quarks 
(assuming Sb ~ 0.7). This ultimately yields something more than 44 events per year. In 
addition, one should expect most of the signal events to lie in the detector acceptance region, 
since leptons and jets originate from decays of heavy objects. 

The same signature could well be exploited in the case of the iii2h final state. Here, the 
production cross section is smaller than for iiilh production, about 40 fb, so is the ^2 xfb 
decay rate as compared to the ti — > xfb one (see Tab. However, the final number of 
events per year is still quite large: about 14 after having already multiplied by e^. 

As for the kinematics of these two signatures, we may remark that they have peculiar 
features that should help in their selection: a not too large hadronic multiplicity, six jets in 
total, each rather energetic (in fact, note that m^± — m^o ^ 58 GeV and m^^ ~ 378 GeV ^ 
~ 273 GeV ^ m^± ^114 GeV), so that their reconstruction from the detected tracks 
should be reasonably accurate; high transverse momentum and isolated leptons to be used 
as trigger; large -Emiss to reduce non-SUSY processes; four tagged 6-jets that can be exploited 
to suppress the + light jet' background from QCD, and one bb pair resonating at the 
h mass, m,h ~ 90 GeV. Moreover, the 'irreducible' background from titlZ events has been 



shown in Ref. [jT9| to be under control, even when rrih ~ Mz, as it is the case here. 
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Particle 
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ii 


76% 


xtb 




19% 


x\t 




57% 


xtb 




24% 


xtb 




11% 


xlt 


h 


90% 


bb 




5% 


1+1- 




0.0003% 


77 



Particle BR Decay 



t 



xt 



xl 



33% 
3_3% 
11% 
11% 
11% 



3_0% 
3_0% 
14% 
14% 
14% 



28% 
14% 
14% 
14% 



udb 
csb 
e^vb 
fi^ub 
T^ub 



Xics^ 

x\r^y^ 



Table 5: Dominant decay channels and BRs of the final state (s)particles in gg — >■ QxQph, 
q = t and x, x' = 1, 2, for Mq = M1/2 = 130 GeV, Aq = 0, tan/3 = 3 and sign(/i) < [0. 
Via ofF-shell W+. 
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Particle BR Decay 



t2 



hi 



h 



H 



A 



94% 



40% 
26% 
16% 
7% 



6_1% 
32% 



4_2% 
31% 
18% 



94% 
6% 



94% 
6% 



94% 
6% 



9_1% 

5% 



Xih 



xtb 
xtb 
biW+ 



X\h 



Xlh 
Xlb 
X'2b 



bb 



T ' T 



bb 



T ' T 



bb 



T ' T 



X\Xl 



Particle BR Decay 



Xi 



X2 



xl 



xl 



xl 



100% 



24% 
15% 
11% 
10% 



100% 



24% 
24% 
9% 



2_3% 
2_3% 
6% 
6% 



100% 



xlw+ 
xtz 

xtA 



xtw- 

XiW+ 

xlz 



xtW' 

XiW^ 
X\h 

xlz 



Table 6: Dominant decay channels and BRs of the final state (s)particles in gg — >■ q^q*,X, 



X = h,H, A, H^, q = t,b and x, x' 
and sign(/i) < [[T0 |. 



1,2, for Mo = M_ 



1/2 



130 GeV, Aq = 0, tan/3 = 40 
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4.1.6.2 Large tan/5 regime 

In the large tan (3 regime there is a variety of cross sections which can be significant: titJX, 
iii*X, biblX, hhlX, iii*A, bib* A, and iiblH+ , where X = h,H. For reasons of 

space, however, we only focus our attention to one signature for the Higgs states not yet 
considered, the one arising from the dominant production channel in all cases, with the 
only exception of the charged Higgs bosons. In fact, we will neglect analysing here their 
decay patterns, as we have already mentioned the poor effectiveness of the charged Higgs 
production modes both as discovery channels and probes of the underlying M-SUGRA model. 

In the case of heavy scalar Higgs bosons, we consider the final state bib\H. This yields 
a cross section of 14 fb for tan (3 = 40. A possible decay chain is the one below. 

h hi H 

J' J' J' 

Xi + b X? + 6 b + b 

That is, a rather simple final state made up by four 6-quarks and missing energy. The BR 
of this sequence is about 9% (see Tab. ^. Therefore, at high luminosity, one obtains 129 
events per year, before heavy flavour identification. 

The main background is certainly ordinary QCD production of four jets. However, the 
requirement of tagging four 6-jets would reduce the latter considerably. Furthermore, if the 
mass of the heavy scalar Higgs boson is known, then one could impose that two 6-quarks 
reproduce tuh ~ 121 GeV within a few GeV (say, 5) in invariant mass. Finally, given the 
enormous mass difference — rub ~ 250 GeV, compared to the rest mass of the LSP, 
m^o f» 51 GeV, one should expect, on the one hand, a large amount of missing energy, and, 
on the other hand, all four 6-jets to be quite hard, both aspects further helping to reduce 
the QCD noise. In the end, some 32 events could well be detected annually, having already 
accounted for the overall 6-tagging efficiency ef = 0.25. 

For the case of the pseudoscalar Higgs particle, we choose the final state titg^, which 
has a cross section of about 79 fb at tan/3 = 40. A possible signature could be^: 









A 


i 


i 




i 




Xi 


+ b 


b + b 


i 


i 










+ V 




i 


i 


+ xS 











Here, the final state is made up by four 6-quarks and two r's, plus missing energy as usual. 
The decay fraction is 23% (again, see Tab. |^). That is, 1814 events per year before tagging 
6's and r's. 

^Additional examples can be found in Ref . [16| . 
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The most dangerous backgrounds are probably Z + 4 jet production and ttbb events. The 
first can be rejected by asking, e.g., M^+^- 7^ Mz, if r's are reconstructed. In addition, 
both background processes have (at least) two 6-quarks quite soft. As for the signal, given 
that the lightest chargino mass is much smaller than the stop ones (in fact, m-± ^ 93 GeV 
whereas m^^ ^ 248 GeV and m^^ ~ 388 GeV), all 6's are naturally energetic and two of them 
also peak at rriA ~ 120 GeV. Thus, to require all Mbb invariant masses sufficiently large with 
one close to the A mass should help in enhancing considerably the signal-to-background 
rates. Requiring large missing energy would help further, especially against tibb events. 
More difficult is to discern differences in the r behaviours (though, notice that rrif^ ^ 76 
GeV ^ rrir). For ef = 0.25, and assuming leptonic decays of both r"*" and r~ into electrons 
and/or muons, one finally gets something of the order of 110 signal events per year. 

4.2 Heavy mass spectrum 

An attempt to summarise our findings is made in Fig. ^, where the most relevant cross 
sections (see upper frame) for the light mass regime (see lower frame) are plotted for a 
choice of Mq, M1/2, Aq, tan/3 and sign(/i) which refiects their hierarchal order seen over 
most of the M-SUGRA parameter space discussed so far. 

However, to assume that only light squark and Higgs masses can induce sizable cross 
sections in events of the type (|l]) would be wrong. This is a sufficient condition for many 
channels, but not a necessary one. For example, even for very large universal masses, one 
could find a value of the soft trilinear coupling small enough to overcome the loss of signal due 
to propagator and phase space effects. In fact, as repeatedly shown in the previous Section, 
most of the cross sections considered here grow quickly as Aq becomes negative. Fig. § makes 
eloquently this point (top insert), for the choice Aq = —900 GeV, well consistent with the 
bounds imposed by the charge and colour breaking minima. There, we have adopted a very 
large value for Mq, i.e., 500 GeV, and varied M1/2 between 100 and 500 GeV. The choice 
of a large tan/? value, i.e., 35, is necessary to obtain detectable rates, except in those cases 
involving ti and h in the same event. In contrast, that of a negative sign(/i) never is. The 
squark and Higgs masses produced by the above combinations of M-SUGRA parameters can 
be found in the bottom frame of Fig. ^. 

There are only a few production channels which survive the strong phase space suppres- 
sion arising in the heavy mass regime and yield cross sections above 1 fb. Among these, 
other than those already encountered titlh, tit^^h, bib^h, one notices the appearance of chan- 
nels which had negligible rates in the low mass regime, notably 6162^ (compare to Fig. 
In this specific instance, the effect is due to the enhancement induced by the onset of the 
62 — ^ bih decay mode in gg — 62^2 events. Even the other two combinations tit^H and 
tifog-f^", which had a rather low profile in the light mass regime, now compete more closely 
with the dominant modes. Finally, notice that pseudoscalar Higgs boson production is no 
longer significant in this mass regime, in line with the results presented in Ref. [p^] . 

In practise, if Mq = 500 GeV, events involving light CP-even Higgs bosons could be 
detected up to M1/2 = 400 GeV or so, in either mode tii\h or fit^^/i- Final states of the type 
bib^h have surprisingly large rates if M1/2 is below 220 GeV. The maximum reach in M1/2 
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via biblh, tih*2H and tit\H is instead 220, 180 and 140 GeV, respectively. 

This is just one example where a new phenomenology of squark-squark-Higgs events 
arises for rather heavy Mq and M1/2 masses. We have found several such combinations, and 
linked them to the fact that Aq ought to be significantly large and negative, tan [3 close to 
nit /nib, but with small dependence on sign(/i). 



5 Conclusions 

In summary, we have studied neutral and charged Higgs boson production in association 
with all possible combinations of stop and sbottom squarks at the LHC, in the context of 
the SUGRA inspired MSSM. Our interest in such reactions was driven not only by the fact 
that they can act as new sources of Higgs particles but also because they carry a strong 
dependence on the five inputs of the SUSY model, so that they can possibly be used to 
constrain the latter. In a sense, this note (along with Ref. |16|) completes previous analyses 
on the subject ^7\, |TB|, |T^, where the emphasis was mainly put on the usefulness of the above 



kind of reactions as Higgs production modes and the attention consequently restricted to 
the case of light squark and Higgs masses. 

We have found that the cross sections of many of these processes should be detectable at 
high collider luminosity for not too small values of tan j3. Indeed, their production rates are 
strongly sensitive to the ratio of the VEVs of the Higgs fields, this possibly allowing one to 
put potent constraints on such a crucial parameter of the MSSM Higgs sector. Furthermore, 
also the trilinear coupling Aq intervenes in these events, in such a way that visible rates 
would mainly be possible if this other fundamental M-SUGRA input is negative. (Indeed, to 
know the actual value of tan f3 from other sources would further help to assess the magnitude 
of Aq.) As for the sign of the Higgsino mass term, sign(yu), it affects the phenomenology 
of such events in one or two channels only, so that it can easily evades the imposition of 
experimental bounds. Finally, concerning the remaining two parameters (apart from mixing 
effects) of the M-SUGRA scenario, one must say that Mq need not be small (it could be 
as large as 500 GeV) and that M1/2 is enough to be below 220 GeV in order to guarantee 
sizable cross sections in many cases. 

In a few representative examples, we have further investigated the decay phenomenology 
of these reactions, by discussing some possible signatures, their rates (of the order of tens to 
hundreds of events per year at high luminosity) and peculiar kinematics, as opposed to the 
yield of ordinary, non-SUSY backgrounds. 

In conclusion, we believe these processes to be potentially very helpful in putting drastic 
limits on several M-SUGRA parameters and we thus recommend that their phenomenology is 
further investigated in the context of dedicated experimental simulations, which were clearly 
beyond the scope of this note. In this spirit, we have derived compact analytical formulae 
of the relevant production MEs, that can easily be incorporated in existing MC programs. 
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Appendix A 

In this additional Section we present in analytic form the MEs adopted in calculating all our 
processes. We identify the external particles as follows: 

9{Pi. Ai) + g{p2, A2) + '?x'(P4) + $(P5), (A.l) 

where Pi, i = 1, . . . 5 are the four-moment aQ and Xj, j = 1, 2, are the helicities of the QCD 
vector bosons. 

It is convenient to rearrange the amplitudes corresponding to the graphs in Fig |l] in 
terms of their colour structure. For example, it is trivial to recognise the existence of only 
two combinations of colour matrices, namely, t^ctcb = (t^t^)ab and t^c^cb = (t^t'^)ab, where 
A,B = 1, ... 8 are the gluon and a,b,c = 1, ... 3 the squark colour indices. This is imme- 
diate for graphs 1 to 6, as one can realise by explicitly writing down the QCD Feynman 
rules. In addition, the triple-gluon diagrams, graphs 7 and 8, are proportional to the anti- 
commutator [^^,1^], whereas those involving quartic couplings, graphs 9 and 10, depend on 
the commutator {t^,t^}. 

Therefore, the total amplitude of processes of the type ( |A.1| ) can be written as 

Ai%,b = {t^tXbT}'^ + {t^tX,Ti'\ (A.2) 

The two subamplitudes t/'^^, i = 1,2, are obtained as follows 

3 8 10 



^ 5^M,^^> + 5^Mf>+5^Mf>, (A.3) 

i=l i=7 i=9 

6 8 10 

i=4 i=7 1=9 

where the m}^\ i = 1, . . . 10, are the original Feynman amplitudes associated to the graphs 
in Fig. 0, but deprived of their colour structure (and couplings, see eq. ( [A.4| ) below)|^. 

This way, the total amplitude squared, summed/averaged over the final/initial spin and 
colours, can be expressed in terms of only two colour factors, as 

\M?(99 - ^.$*) = I W ^tt T^''Tj>^'C,,, (A.4) 

{A}=± i=l j=l 

with g"^ = AtcOs, g^r = ATcaem/siy and where Cij is a (2 x 2) colour matrix with elements 

Cu s C22 = i(^-2iVc + iVil) = y, (A.5) 



"'^^For our purposes, we take the initial state momenta, pi and p2, as incoming and the final state ones, ps, 
P4 and P5, as outgoing. 

^^The notation {A} refers cumulatively to the helicities of the two incoming gluons, i.e., with i — 1,2. 
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as usual being Nc = 3 the number of colours in QCD. In eq. (|A.4|) , \q,q^q' ^ represents the 
strength of the squark-squark-Higgs vertex involved, as described in Sect. ^ (apart from an 
overall phase and the factor gw)- 



The ten amplitudes m/^^ are simply 



Ml 



{A} 



M 

M. 

M 



{A} 
2 

{A} 

3 

{A} 
4 

{A} 
5 

{A} 

{A} 
7 

{A} 

3 

{A} 
3 

{A} 
10 



4£i(Ai 
4£i(Ai 
4£i(Ai 

4£i(Ai 
4£i(Ai 

£^12(Al, A2 



Pi e2{\2) ■ {P3+P5)/Pu/P35, 
iP4+P5) £^2(A2) ■P3/P23/P45, 

Pa £2(A2) ■P3/-P14/-P23, 

(P3 +P5) £^2(A2) -^4/^4/ As, 
Pa ^r2(A2) • (P4 +P5)M3/i^45, 
P3 £^2(A2) -^4/^4/^3, 
^ ^-2p4 +Pl +P2)/i^35, 



£:i2(Al, A2) • (+2p3 -pi 

^l(Al)-£2(A2)/P35, 

£l(Al)-52(A2)/P45, 



P2)/P. 



45, 



(A.6) 



where we have introduced the propagator functions 



Pl4 = 


(Pl 


-P.? 


-Ml, 


-P24 = 


{P2 


-P.? 




-P23 = 


{P2 


-P3? 




Pl3 = 


{Pl 


-P3? 




P35 = 


{P3 


+p,r 


-M|,+iM,,r,, 


-P45 = 


{Pi 


+ P.f 





(A.7) 



the gluon polarisation vectors |^T|, i = 1, 2 



±) = ^[Ter(A. = l)-i5f(A, = 2)], 

^r(A. = i) = mpJr\^.P-PiP\pi-pf). 
.f(A. = 2) = (pn-^(o,-pr,pr,o). 



(A.8) 



with 



T 

Pi 



+pf , 
pf +pf +pf, 



(A.9) 
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and their contraction over the triple-gluon vertex times the gluon propagator 

^^2(Ai,A2) = ^{(pi-P2)%i(Ai)-£2(A2) (A.IO) 
-n.2 

+ [(P2 +P12) • £l(Al)]£^(A2) - [(Pl +P12) • £2(A2)]£^(Ai)}, 

with P12 = = (Pi +P2)^- 
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graph 9 




graph 10 



Figure 1: Feynman diagrams contributing at lowest order to processes ( |A.1|) . An helical line 
refers to a gluon whereas a dashed one symbolises both a squark and a Higgs boson. 




Figure 2: Resulting masses from for the third generation of squarks versus tan (3 for 
three values of Aq and for both positive (dot-dashed) and negative (solid) fi. 




Figure 3: Resulting masses from for the five Higgs bosons versus tan (3 for three values 
of Aq and for both positive (dot-dashed) and negative (solid) fi. 
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Figure 4: Total cross section of gg qi,2Qi,2h processes (with q = t,b) as a function of 
tan/3 for the characteristic input values Mq = 150 GeV and M1/2 = 150 GeV. Both positive 
(dot-dashed) and negative (solid) /i as well several Aq contour lines are shown. The symbol 
"x" is used to indicate parameter areas forbidden by direct Higgs boson searches. 



H 



bi b2 H 



b 



10 



10 



10 



/j,<0 GeV (solid) 
/x>0 GeV (dotdash) 
M^/2=150 
Mq=150 



10" 



10" 




.J2 

b 



IQJ 



10^ 



10" 



10" 



10" 



/i<0 GeV (solid) 
/x>0 GeV (dotdash) 
M^/2=150 GeV 
Mq=150 GeV 



Vs=14 TeVi 




b2 b2 H 



H 



J3 

b 



b 



10^ 



10^ 



10" 



10" 



10" 



= 1 1 1 1 1 1 1 1 1 1 
: /j,<0 GeV (solid) 


1 1 1 1 z 

Vs=14 TeV: 


IJ>0 GeV (dotdash) 








: Mq=150 GeV ^ 




jl /\ \ 1 1 1 1 1 1 





b 



10 20 

tan;8 



30 




10"= 



IQJ 



/j,<0 GeV (solid) 
/j,>0 GeV (dotdash) 
M^/2 = 1S0 "^eV 
Mq=150 GeV 



Vs=14 TeVi 




tan/S 



^2 t2 H 



l.OO 
0.50 

0.20 
0.10 
0.05 

0.02 
0.01 



^t<0 GeV (solid) 
,((-^0 GeV (dotdash) 



= 150 GeV 



Ma=150 GeV 



Vs=14 TeV- 




tan/S 



10 20 

tan/S 



30 



Figure 5: Total cross section of gg qi,2(li,2H processes (with g = t, 6) as a function of 
tan/3 for the characteristic input values Mq = 150 GeV and M1/2 = 150 GeV. Both positive 
(dot-dashed) and negative (solid) fi as well several Aq contour lines are shown. The symbol 
"x" is used to indicate parameter areas forbidden by direct Higgs boson searches. 
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Figure 6: Total cross section of gg qiq2A processes (with g = t, 6) as a function of tan /? 
for the characteristic input values Mq = 150 GeV and M1/2 = 150 GeV. Both positive (dot- 
dashed) and negative (sohd) n as well several Aq contour hues are shown. The symbol "x" 
is used to indicate parameter areas forbidden by direct Higgs boson searches. 
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Figure 7: Total cross section of gg processes (with gC) = t, h) function of 

tan/5 for the characteristic input values Mq = 150 GeV and M1/2 = 150 GeV. Both positive 
(dot-dashed) and negative (solid) ix as well several Aq contour lines are shown. The symbol 
"x" is used to indicate parameter areas forbidden by direct Higgs boson searches. 




Figure 8: The most significant cross sections which survive the choice of a hght SUSY 
spectrum (above), alongside the values for the masses entering the corresponding production 
processes (below). Shaded regions indicate areas excluded by direct searches. 




Figure 9: The most significant cross sections which survive the choice of a heavy SUSY 
spectrum (above), alongside the values for the masses entering the corresponding production 
processes (below). Shaded regions indicate areas excluded by direct searches. 



